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■ Abstract. 

We model the responses of Tissue-Equivalent Proportional Counters (TEPC) 
to radiation fields of therapeutic 12 C beams in a water phantom and to quasi- 
monoenergetic neutrons in a PMMA phantom. Simulations are performed with the 
Monte Carlo for Heavy Ion Therapy (MCHIT) model based on the Geant4 toolkit. The 
shapes of the calculated lineal energy spectra agree well with measurements in both 
cases. Measurements on the axis of a narrow 12 C beam with its width smaller than 
the TEPC diameter are studied in detail. The relation between LET and measured 
frequency- mean lineal energy £// is discussed. It is found that the choice of the nuclear 
fragmentation model used in MCHIT simulations has a relatively small influence on 
the calculated total lineal energy spectra. However, the shapes of the spectra measured 
with 12 C beams are better described by MCHIT calculations which take into account 
the production of 5-electrons. The validation of MCHIT with neutron beams gives us 
confidence in estimating the contributions to lineal energy spectra due to secondary 
neutrons produced in water by 12 C nuclei. As found, the neutron contribution at 10 cm 
distance from the beam axis amounts to ~ 40% close the entrance to the phantom and 
decreases to ~ 25% at the depth of the Bragg peak and beyond it. The presented 
results can help in evaluating biological out-of-field doses in carbon-ion therapy. 

Keywords: theory and algorithms in physics of heavy-ion therapy, microdosimetry, 
Monte Carlo applications, simulation 
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1. Introduction 

At present time the heavy-ion therapy of cancer is one of the most advanced 
methods in radiation therapy (Schardt et al. 2010, Elsaesser et al. 2010, Durante 
& Loeffler 2010). An elevated biological effectiveness of accelerated nuclei provides 
advantages in treatment of radioresistant solid tumors, but also requires a thorough 
treatment planning to reduce undesirable impact of radiation on healthy tissues. This 
should include, in particular, a realistic description of the production of secondary 
nuclear fragments. Such secondary particles deliver dose to the tissues located farther 
than the Bragg peak for primary nuclei and also around the primary beams. 

Recently there were several experimental (Haettner et al. 2006) and theoreti- 
cal (Gudowska et al. 2004, Bohlen et al. 2010, Pshenichnov et al. 2010) studies which 
confirmed the need to consider nuclear fragmentation reactions in heavy-ion therapy. In 
particular, as shown by measurements with a A00A MeV 12 C beam, about 70% of beam 
nuclei undergo fragmentation reactions (Haettner et al. 2006). Apart from the effect 
of primary beam attenuation, this means that the dose around the beam is delivered 
by light charged fragments, (e.g. protons and helium nuclei) and also by neutrons. In 
particular, neutrons may propagate large distances from their production points before 
they initiate secondary nuclear reactions. 

As demonstrated by studies with various transport codes, such as SHIELD- 
HIT (Gudowska et al. 2004), FLUKA and Geant4 (Bohlen et al. 2010, Pshenichnov 
et al. 2010), the measured attenuation of the primary carbon beam and production of 
heavier fragments are generally well reproduced by the theory. However, as follows from 
the same studies, there are problems with quantitative description of the yields, angular 
distributions and energy spectra for all kinds of fragments produced in fragmentation 
reactions in water. In particular, the yields of alpha-particles are underestimated. This 
may be related to neglecting the cluster structure of light nuclei, which leads to the 
underestimation of decays of excited nuclei into channels containing alphas. It is unlikely 
that such deficiency will be eliminated soon, as this will require combining dynamical 
and statistical models of nuclear reactions with sophisticated nuclear structure models 
into a unified, but still computationally effective tool. 

Nuclear fragmentation models calculate the yields of secondary fragments as 
functions of their charges and velocities. However, the impact of specific particle 
species on living cells is rather correlated with their linear energy transfer (LET) 
values. The LETs of particles with different charges Z, but similar velocities in first 
approximation are proportional to their Z 2 (Taddei et al. 2006). At the same time 
ions of different velocities and charges may have similar LET (Guetersloh et al. 2004). 
Therefore, the reliability of nuclear fragmentation models, in particular used in Geant4 
calculations (Bohlen et al. 2010, Pshenichnov et al. 2010), should be also evaluated with 
respect to calculations of LET or related quantities. 

The biological action of radiation is commonly characterized by microdosimetric 
quantities (International Commission of Radiation Units and Measurements (ICRU) 
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1983, Rossi 1990, Lindborg & Nikjoo 2011) such as lineal energy y, its probability 
density f(y) and dose probability density d(y). Frequency- mean, y~f, and dose- mean, 
yd, lineal energies are defined as the first moments of the corresponding distributions, 
see section 2 for their definitions. These distributions are routinely measured by means 
of Tissue Equivalent Proportional Counters (TEPC) (Waker 1995, Waker 2006). 

In a complex radiation field different particles can eventually contribute with similar 
y. Therefore, the resulting y-distribution is build as the sum of contributions from 
primary and all secondary particles passing TEPC. Measurements and calculations 
of microdosimetric quantities related to heavy-ion therapy were reported by several 
authors. In particular, microdosimetric quantities were used to characterize the 
biological dose from therapeutic 12 C beams (Kase et al. 2006, Kase et al. 2011, Martino 
et al. 2010). Similar measurements were performed for the first time for 7 Li 
beams (Martino et al. 2010). Simulations with the PHITS code were done for a wall- 
less TEPC (Tsuda et al. 2010). The energy deposition inside and around 7 Li and 
12 C beams in water were calculated with SHIELD-HIT and FLUKA codes (Hultqvist 
et al. 2010, Taleei et al. 2011), but without implementing the exact geometry of the 
TEPC detector and, respectively, without calculating ^-distributions. The FLUKA 
code was also used to simulate responses of TEPC to various ions (Bohlen et al. 2011), 
including ^-distributions for TEPCs located on the beam axis inside a water phantom 
irradiated by 12 C beam (Bohlen et al. 2012). 

In our previous papers (Pshenichnov et al. 2005, Pshenichnov et al. 2006, 
Pshenichnov et al. 2007, Pshenichnov et al. 2008, Pshenichnov et al. 2010, Mishustin 
et al. 2010) we validated the Geant4 toolkit (Agostinelli et al. 2003, Allison et al. 2006) 
for use in heavy-ion therapy simulations. The depth-dose distributions (Pshenichnov 
et al. 2006) and production of positron-emitting nuclei (Pshenichnov et al. 2007, 
Pshenichnov et al. 2008) were calculated by means of the Monte Carlo for Heavy Ion 
Therapy (MCHIT) model based on the Geant4 toolkit. Special effort was made to 
evaluate the quality of nuclear fragmentation models (Pshenichnov et al. 2010, Mishustin 
et al. 2010) by comparing theoretical predictions with experimental data on fragment 
yields (Haettner et al. 2006). 

As mentioned in a renowned publication on microdosimetry (International 
Commission of Radiation Units and Measurements (ICRU) 1983), page 38: 
"Experiments and calculations of microdosimetric spectra play complementary roles. 
The degree of agreement between experiment and theory serves as a test of the validity 
of both." Following this approach we extend the MCHIT model to calculations of 
microdosimetric quantities suitable for estimating the biological effects of accelerated 
nucleons and nuclei. The results of the calculations are compared with recent 
data (Martino et al. 2010) obtained with a therapeutic 300A MeV 12 C beam and with 
other experimental data. To the best of our knowledge the present study is the first one 
which uses the Geant4 toolkit for calculating microdosimetric quantities by a detailed 
modeling of TEPCs located in a phantom both inside and outside of a therapeutic 
12 C beam. This work complements the above-mentioned theoretical and experimental 
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studies by considering both in-field and out-of-field ^-distributions in a consistent way 
and especially by evaluating the contribution of secondary neutrons to out-of-field doses. 
Our computational methods can be useful to assess the biological effects of secondary 
neutrons produced in carbon-ion therapy (Newhauser & Durante 2011). 

In order to study the conditions of microdosimetry measurements performed at 
the GSI Helmholtz Center for Heavy Ion Research with therapeutic pencil-like 12 C 
beams (Martino et al. 2010) we also calculated the number of tracks and total track 
length for primary and secondary particles inside a TEPC irradiated by such a focused 
beam. This helps to determine £// and y d in the cases when a TEPC is not randomly 
traversed by particles and (or) several particles traverse the TEPC in a single event. 

In section 2 we begin our study by presenting the basic principles of microdosimetry 
measurements, definitions of microdosimetric variables and distributions. Section 3 
briefly describes the Geant4 physics models used in MCHIT. In section 4 we describe the 
modeling of experimental set-up and specific computational techniques used to speed-up 
Monte Carlo calculations of microdosimetric quantities with MCHIT. Our calculations 
of microdosimetric distributions for neutron beams in polymethilmethacrylate (PMMA) 
are presented in section 5 in order to prove the validity of MCHIT for estimating the 
contributions of secondary neutrons produced by 12 C beams in section 7. The central 
part of our study is presented in section 6, where the topology of energy deposition to 
TEPCs placed in a water phantom irradiated by a pencil-like monoenergetic 12 C beam is 
investigated in detail. Calculations of microdosimetric distributions, their corresponding 
moments £//, yd and total dose for the 12 C beam are also presented in section 6. This 
section further contains the discussion of relation between y~f and LET for TEPCs placed 
either in water or in air and irradiated by 12 C beams of various profiles and energies. 
The results of our study are summarized in section 8. 

2. Basics of microdosimetry method and relevant measurements 

It is known that the impact of charged particles on eukaryotic cells is defined not only by 
the energy imparted to the cell volume (i.e. by the absorbed dose), but also by the total 
number and spatial correlations of local energy deposits to DNA molecules inside the cell 
nucleus. Such deposits lead to spatially distributed DNA lesions caused mainly by free 
electrons which surround a track of a charged particle. The pattern of stochastic impacts 
on DNA resulting in single, double or more complex breaks of DNA strands essentially 
depends on the projectile charge, mass and velocity. It is defined by the number of 
electrons per unit of the track length as well as on their energies. All these characteristics 
differ for photon, electron, proton and nuclear beams. In particular, carbon nuclei 
successfully used for particle therapy of cancer are characterized by elevated biological 
effectiveness in the Bragg peak region as compared to photons and protons (Elsaesser 
et al. 2010, Durante & Loeffler 2010). Thus, these biological effects of ion beams should 
be taken into account while developing the treatment planning systems used in particle 
therapy (Kase et al. 2006, Kraemer & Durante 2010, Kase et al. 2011). 
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Since DNA molecules are typically confined inside cell nuclei, which are structures 
of a few-micrometer size, it is desirable to study energy deposition patterns at the 
micrometer scale. In particular, the in-beam microscopy is used in experiments on 
irradiation of living cells. With this tool, the tracks of charged particles can be well 
explored, if visualized by biological markers for DNA damage (Tobias et al. 2010). 

Since many years measurements with Tissue-Equivalent Proportional Counters 
(TEPC) (International Commission of Radiation Units and Measurements (ICRU) 1983) 
are considered as a practical tool for studying patterns of energy deposition to 
micrometer-size objects by various kinds of radiation. The TEPC considered in this 
paper has a chamber with an inner diameter of 12.7 mm filled with low-pressure 
tissue-equivalent (TE) gas (Martino et al. 2010). It emulates a tissue-equivalent 
region of a micrometer size. In this way the energy deposited to the sensitive gas 
volume of this device can be related to the energy that would be deposited to a cell 
nucleus (International Commission of Radiation Units and Measurements (ICRU) 1983). 
As demonstrated by several authors, see e.g. (Rossi 1990, Lindborg & Nikjoo 2011), 
microdosimetric variables are directly related to the radiation quality factor used to 
quantify risks from various kinds of radiation. 

Due to the stochastic nature of particle transport in media, the amount of energy 
e delivered to a sensitive volume representing a cell nucleus changes from one event 
to another. The value of e in each energy deposition event can be measured by a 
TEPC (Waker 1995, Waker 2006, Lindborg & Nikjoo 2011), as the measured electronic 
signal is expected to be proportional to e. This technique allows one to obtain the 
probability distributions for lineal energy y = e/l, where I is the mean chord length of 
the sensitive volume used as a detector. If one assumes that particles representing a 
homogeneous radiation field traverse a spherical detector randomly at various impact 
parameters, then I — 2 /3d. The specific energy z = e/m can be introduced as 
another stochastic quantity, where m is defined as the mass of matter in the sensitive 
volume (International Commission of Radiation Units and Measurements (ICRU) 1983). 

There were several experimental (Guetersloh et al. 2004, Taddei et al. 2006) and 
theoretical (Nikjoo et al. 2002, Taddei et al. 2008) studies of TEPC responses to heavy- 
ions. In these measurements and calculations the dependence of e on the impact 
parameter of the ion track was investigated in detail. As demonstrated by Taddei and co- 
authors (Taddei et al. 2008), calculations based on the Geant4 toolkit can successfully 
reproduce the performance of a TEPC (Far West Technology Inc., model LET-1/2) 
irradiated by 4 He, 12 C, 16 0, 28 Si and 56 Fe nuclei. 

Since lineal energy y varies from one event to another, the frequency- 
mean lineal energy yf, which characterizes the probability density f{y), can be 
introduced (International Commission of Radiation Units and Measurements (ICRU) 



The value of yf is the first moment of f(y) and serves as a measurable approximation 



1983): 
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for LET. In this work, following the conventions accepted by Martino et al. (Martino 
et al. 2010), we consider f(y) as the probability density of y per beam particle (carbon 
ion). This differs from the usual definition of f(y) per energy deposition event, see 
e.g. (International Commission of Radiation Units and Measurements (ICRU) 1983), as 
accepted by Taddei et al. (Taddei et al. 2008) (section 6.4) and also in measurements 
with neutron beams in PMMA (section 5). As explained below, yj can be used to 
calculate the dose, see Eqs. (3) and (4). In the following also ^(^-distributions, with 
d(y) = yf(y)/y~f defined as the dose probability density, will be considered. The dose- 
mean lineal energy y d is calculated as the second moment of the /(^-distribution divided 
by its first moment: 

1 .2 . 



y d = / yd(y)dy = - / y 2 f(y)dy . (2) 
Jo Vf Jo 

The value of yd is related to the quality of radiation. 

The mean specific energy z can be calculated from microdosimetric variables for a 

spherical sensitive volume with the diameter d and density p as: 

(3) 

m m 7vd z p 

In the case the volume is traversed by a single particle, the mean specific energy z is equal 
to the total dose D, and hence can be calculated from Eq. (3). For a sensitive volume 
filled with gas the total dose D is given by a useful formula (International Commission 
of Radiation Units and Measurements (ICRU) 1983) resulting from Eq. (3): 
^ 0.204 

D = — V, ■ (4) 

Here d corresponds to the diameter of the sensitive volume expressed in /iin, y/ in 
keV//zm and D is obtained in Gy. 

Due to its compactness the TEPC device can be located at various positions 
inside and outside of phantoms irradiated by therapeutic beams. The microdosimetry 
measurements with a monoenergetic 400A MeV 12 C beam (Endo et al. 2005) were 
made together with the identification of charges of secondary fragments by scintillation 
counters by the time-of-flight (TOF) method. In this way the above-described 
microdosimetry technique supplements the measurements of fragment yields by particle 
identification made outside the phantom with more bulky detectors, see e.g. (Haettner 
et al. 2006). 

Further measurements were performed with phantoms irradiated by 290A MeV 
12 C beam with a traditional walled TEPC (Endo, Takada, Onizuka, Tanaka, Maeda, 
Ishikawa, Miyahara, Hayabuchi, Shizuma & Hoshi 2007), also with fragment charge 
identification, and later (Tsuda et al. 2010) with a wall-less counter without tagging 
fragments. Recently the results of microdosimetry measurements with 290A MeV 
Spread Out Bragg Peak (SOBP) carbon-ion beam were also reported (Endo et al. 2010, 
Kase et al. 2011). In all these measurements the TEPCs were located at the beam axis. 

Microdosimetry measurements on the beam axis as well as off-axis in a water 
phantom irradiated by a pencil-like 12 C beam were performed at GSI (Martino 
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et al. 2010). Similar studies at carbon-ion and proton radiotherapy facilities with passive 
beam delivery were also performed at TEPC positions outside the treatment field (Yonai 
et al. 2010). The collected microdosimetry data challenge the theoretical methods aimed 
to describe particle transport in tissue-like media. 

3. Physics models used in MCHIT simulations 

For the purposes of the present study we further developed our Monte Carlo model 
for Heavy-Ion Therapy (MCHIT) (Pshenichnov et al. 2006, Pshenichnov et al. 2007, 
Pshenichnov et al. 2008, Pshenichnov et al. 2010, Mishustin et al. 2010) to simulate 
microdosimetry measurements with TEPC devices. The model is currently based on 
the Geant4 toolkit (Agostinelli et al. 2003, Allison et al. 2006) of version 9.4 with 
patch 02. Input parameters for a MCHIT run are provided via a set of user interface 
commands which define the beam particle, the energy spread of the beam, its spot size, 
angular divergence, as well as the dimensions and composition of the phantom used in 
a particular experimental setup. 

A detailed description of the physics models included in the Geant4 toolkit is given 
in the Geant4 Physics Reference Manual (Geant4 Physics Reference Manual 2011). A set 
of models which are relevant to a particular simulation problem should be activated by 
the application developer. Such models are usually grouped in a physics list. In order to 
involve certain models into simulation, one can either use the so-called predefined physics 
lists, or implement customized physics lists, or even use a combination of these two 
options. The predefined physics lists are provided by Geant4 developers and distributed 
along with the Geant4 source code. It is convenient to use separate physics lists for 
electromagnetic and hadronic physics. Besides the user can activate or de-activate a 
specific physical process, the low-energy neutron capture, for example, in a MCHIT 
simulation via corresponding user interface commands. 

3.1. Electromagnetic physics 

The electromagnetic processes are described by means of several predefined physics 
lists of Geant4. The physics list called "Standard Electromagnetic Physics Option 3" 
is recommended by the Geant4 developers for simulations related to particle therapy. 
As demonstrated (Lechner et al. 2010), the measured positions of the Bragg peak for 
carbon nuclei of various energies in water are well reproduced with this physics list. 
The involved physics models simulate the energy loss and straggling of primary and 
secondary charged particles due to interaction with atomic electrons. The multiple 
Coulomb scattering of charged particles on atomic nuclei is also simulated. 

At each simulation step the ionization energy loss of a charged particle is calculated 
according to the Bethe-Bloch formula or interpolated between values listed in a table, 
depending on the particle type and energy (Geant4 Physics Reference Manual 2011). In 
particular, the Bethe-Bloch formula with the shell, density and high-order corrections 
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is applied in G^BetheBlochModel to protons with kinetic energy above 2 MeV. Below 
2 MeV the stopping power parameterizations (International Commission of Radiation 
Units and Measurements (ICRU) 1993) are used for protons in GJ^BraggModel. The 
same methods are applied to alpha particles taking into account the corresponding mass 
scaling with respect to proton. An interpolation of stopping power tables (International 
Commission of Radiation Units and Measurements (ICRU) 2005) is implemented in 
G^IonParametrisedLossModel to calculate electromagnetic energy loss of ions heavier 
than helium with energies relevant to our study, below 1A GeV. The use of another 
predefined physics list, which extends the capability of electromagnetic models to 
produce and transport (^-electrons down to 100 eV, namely G^EmPenelope based on 
the Penelope model (Geant4 Physics Reference Manual 2011), was also considered. 

In order to reduce the CPU time without affecting the accuracy of calculations, 
different cuts for production of electrons were applied in different materials. They are 
listed in table 1. 

Table 1. Cut in range and energy production threshold for electrons in the TEPC 
materials applied in MCHIT simulations with Standard Electromagnetic Physics 
(Option 3) and with the Penelope model 





Standard opt 3 


Penelope 


material 


cut in 


energy 


cut in 


energy 




range 


threshold 


range 


threshold 




(mm) 


(kcV) 


(mm) 


(keV) 


water 


0.1 


85. 


0.1 


85. 


A- 150 TE plastic 


0.01 


17.6 


0.001 


0.3 


TE gas 


9 


0.99 


3 


0.1 



The cut in range for TE gas was defined by the lowest energy of electrons which can 
be transported by the correspond models. For example, in the Standard Electromagnetic 
Physics (Option 3) list the emission of (^-electrons with energies below 990 eV is not 
simulated in TE gas, but their energies are attributed to the local energy deposition 
and finally added to e. This threshold corresponds to 9 mm range of electrons in gas. 
When the Penelope model is used, this cut in range is reduced to 3 mm. 

In a 1.27 mm-thin plastic shell of the TEPC detector the cut in range for electrons 
was reduced to 10 /im. This helps to simulate accurately grazing interactions of particles 
with the plastic shell, as secondary particles produced in plastic can propagate further 
into gas. According to the default settings of Geant4 there were no cuts for production of 
secondary nucleons and nuclear fragments in electromagnetic models, and all produced 
particles (including electrons) were tracked down to zero kinetic energy. 

3.2. Modeling nuclear fragmentation 

A reliable description of the fragmentation of beam nuclei in tissue-like materials is 
crucial for heavy-ion therapy modeling. Nuclear fragments resulting from inelastic 
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nucleus-nucleus collisions modify biological effects of the primary beam. Light nuclear 
fragments are responsible for a tail in the depth-dose distribution which extends 
beyond the Bragg peak of primary nuclei. A customized physics list was implemented 
for the description of hadronic processes. Modifications in this physics list with 
respect to our previous publications concern the modeling of inelastic nucleus-nucleus 
interactions. In the present study we considered two Geant4 models to describe 
the first stage of nucleus-nucleus collisions, namely, the Light Ion Binary Cascade 
(G4BIC) and the Quantum Molecular Dynamics (G4QMD) models. In our earlier 
publication (Pshenichnov et al. 2010) the G4BIC and G4 abrasion models were employed 
to describe such reactions. Recent improvements in the G4QMD code (Koi 2010), 
as well as the conclusions based on the comparison of G4BIC and G4QMD (Bohlen 
et al. 2010), suggest that the G4QMD model can be also successfully used for modeling 
nuclear fragmentation of carbon nuclei. Therefore, in our calculations of microdosimetric 
distributions both options to simulate nuclear fragmentation, G4BIC and G4QMD, were 
used. 

Excited nuclear fragments are frequently produced in addition to free nucleons as 
a result of simulation of a nucleus-nucleus collision event by the G4BIC and G4QMD 
models. Therefore, the G4ExcitationHandler of Geant4 is used to simulate subsequent 
decays of excited nuclear fragments by applying various de-excitation models depending 
on the mass and excitation energy of these fragments. 

The Fermi break-up model (G4FermiBreakUp) is applied to nuclei lighter 
than fluorine. It is designed to describe explosive disintegration of excited light 
nuclei (Bondorf et al. 1995) and it is highly relevant to collisions of light nuclei with 
nuclei of tissue-like materials. As demonstrated (Mishustin et al. 2010) G4BIC linked 
with G4FermiBreakUp better describes the production of lithium and beryllium nuclei, 
as well as secondary neutrons by 12 C nuclei in water and PMMA compared to the option 
of G4BIC linked with the nuclear evaporation model. As demonstrated by Bohlen et 
al. (Bohlen et al. 2010) the combination of G4QMD and G4FermiBreakUp also describes 
data well. Therefore, we used G4FermiBreakUp to de-excite light fragments also in the 
present work. 

For more heavy excited nuclei either evaporation model (Weisskopf & Ewing 
1940) can be used at low excitations (below 3 MeV per nucleon), or the 
statistical multifragmentation model (SMM) (Bondorf et al. 1995) at higher excitation 
energies. The implementation of SMM in Geant4 (G4StatMF) was recently 
validated (Pshenichnov et al. 2010) for studies of nuclear fragmentation reactions in 
extended media. Generally, inclusion of fragmentation reactions helps to describe the 
yield of intermediate-mass fragments. 

3.3. Hits and deposited energy 

We use the primitive scorer classes of Geant4 to calculate physical quantities which 
characterize the TEPC response in each event. The energy e deposited to the TEPC 
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sensitive volume is calculated by means of the G^PSEnergyDeposit scorer. This scorer 
stores a sum of particles' energy deposits to the sensitive volume in each event. 
Employing GJ^PSPassageCellCurrent the number of tracks N that pass through the 
TEPC sensitive volume is calculated. Finally, G^PSPassageTrackLength is used to 
calculate the total track length I inside the TEPC sensitive volume in single and multi- 
particle events. With this scorer the track length is defined as the sum of step lengths 
of the particles inside the volume, and only tracks which traverse the volume are taken 
into account. This means that newly-generated or stopped tracks inside the volume 
are excluded from the calculation of The employed G^PSPassageTrackLength scorer 
is thus insensitive to those secondary electrons which were produced by a fast particle 
inside the TEPC sensitive volume, but stopped inside it. 

4. Modeling experimental set-up 

In this work we calculate microdosimetric distributions and compare them with results 
of two experiments (Nakane & Sakamoto 2001, Martino et al. 2010). In the first 
experiment a PMMA phantom was irradiated by quasi-monoenergetic neutron beams 
at the Japan Atomic Energy Research Institute (JAERI). In the second experiment a 
water phantom was irradiated by carbon and lithium nuclei at GSI. The same type of 
TEPC device (Far West Technology Inc., model LET-1/2) was used by both groups. 
This allowed us to validate our calculational approach for this specific TEPC model first 
with neutron irradiation data (Nakane & Sakamoto 2001), and then extend it to the 
GSI measurements (Martino et al. 2010) with 12 C beam. In this way the contribution of 
secondary neutrons from fragmentation reactions induced by 12 C nuclei can be estimated 
with confidence. 

Experimental results for a TEPC located at 5 cm depth inside the PMMA phantom 
irradiated by 40 MeV and 65 MeV quasi-monoenergetic neutrons were selected from 
JAERI data to validate the MCHIT model. The sensitive region of the TEPC was 
modeled as a sphere filled with a propane-based tissue equivalent gas with 9.03 kPa 
pressure. The device used at JAERI emulated a tissue sphere of 2.07 /iin in diameter. 

In the second experiment (Martino et al. 2010) a water phantom of 30 x 30 x 30 cm 3 
(including side walls made of PMMA) was irradiated by a 300A MeV 12 C beam. In the 
simulations it is assumed that the beam has a concentric spot with a size of 3 mm 
FWHM, an angular spread of 1 mrad both in vertical and horizontal directions, and 
an energy spread of 0.2%. The y<i(t/)-distributions measured at nine points inside the 
water phantom were reported (Martino et al. 2010). They correspond to cm, 2 cm 
and 10 cm radii and to three values of depth: at the plateau of the depth-dose curve, 
at the Bragg peak and in the tail region. In this experiment water was enclosed in 
a container with a 20 mm thick PMMA wall. Carbon nuclei also penetrated through 
a vacuum window made of aluminum, a scintillation detector and a parallel-plate ion 
chamber installed in front of the phantom. These beam-line elements and PMMA 
walls of the phantom were represented in simulations by equivalent water thickness of 
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25.1 mm. This correction was applied to the TEPC positions in order to compare the 
data with simulation results, since calculations were performed with an equivalent wall- 
less water phantom of 30 x 30 x 30 cm 3 . As explained below in section 6.2, a 2 mm 
shift was additionally applied to the TEPC depth values in our calculations with respect 
to experiment in order to describe better the measured microdosimetric spectra at the 
Bragg peak. The TEPC positions used in our simulations are listed in table 2. 

In simulations we have performed for both experiments the actual TEPC dimensions 
and the chemical composition of TEPC materials were implemented as accurate as 
possible. It was assumed that the tissue equivalent gas was propane C 3 H 8 (55%) with 
addition of CO2 (39.6%) and N2 (5.4%). A tissue equivalent plastic shell was defined in 
simulations as a gas container composed of H (10.1%), C (77.6%), N (3.5%), O (5.2%), 
F (1.7%) and Ca (1.9%), in accordance with the properties of the NIST material Shonka 
A-150 tissue-equivalent plastic. The aluminum shell of the TEPC was also introduced, 
although simulation results were found to be rather insensitive to the presence of the Al 
shell because of its small thickness. The central wire located inside the gas sphere was 
not considered in the simulations. Following Taddei et al. (Taddei et al. 2006), it was 
assumed that quite rare events in which a projectile particle interacted with this wire 
are rejected in measurements. We also neglect influence on the simulated tracks due to 
the voltage applied to the internal components of the TEPC. 

Using the axial symmetry of the setup we replaced a single physical TEPC by a 
ring of identical virtual TEPCs located at the same radial distance from the beam axis. 
During each run the number of hits in all counters was scored in a yd(y) histogram 
corresponding to a particular ring. At the end of the run the results were divided by the 
number of virtual counters in the ring and thus became equivalent to a single physical 
counter. With this method one can increase the simulated number of hits for those 
positions of counters where hits are scarce. This is particularly important for TEPCs 
located at large distances, e.g., at 10 cm, from the beam axis. The number of counters 
in each ring and the positions of the rings with respect to the beam axis were chosen 
to avoid any crosstalk impacts on the virtual detectors. This was done to satisfy the 
condition that the particles produced inside one of the virtual counters, or which have 
just traversed it, should not hit other counters. Therefore, calculations for the nine 
TEPC positions listed in table 2 were split into three independent runs. During the run 
"I" three outer rings of TEPCs of 10 cm in radius each consisting of 24 counters were 
placed in the water phantom. In the same run a counter on the beam axis and a ring 
of 2 cm in radius consisting of four TEPCs were placed at 27.71 cm depth. Each of the 
runs "II" and "III" was executed with a ring of 2 cm in radius consisting of four counters 
and with a single counter located at the same depth in the center of the ring. Typically, 
the histories of 10 7 to 5 x 10 7 12 C nuclei traversing the phantom were simulated in each 
run. 

Each counter was modeled as filled with a tissue-equivalent gas with the 
same chemical composition and pressure (12 kPa) as used in the measurements at 
GSI (Martino et al. 2010). The gas was confined in a plastic sphere with a 12.7 mm 
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Table 2. Positions of the TEPC counters inside the water phantom used in MCHIT 
simulations and their labeling in the following text and figures. TEPCs were grouped 
into rings with a number of counters in each ring depending on the ring radius. Each 
TEPC position was modeled in a certain run labeled in the last column. 

radius depth TEPCs position number of run 



(mm) 


(mm) 


notation 


counters 
in the ring 


label 


0. 


52.1 


cm, plateau 


1 


III 


0. 


179.1 


cm, peak 


1 


II 


0. 


277.1 


cm, tail 


1 


I 


20. 


52.1 


2 cm, plateau 


4 


III 


20. 


179.1 


2 cm, peak 


4 


II 


20. 


277.1 


2 cm, tail 


4 


I 


100. 


52.1 


10 cm, plateau 


24 


I 


100. 


179.1 


10 cm, peak 


24 


I 


100. 


277.1 


10 cm, tail 


24 


I 



inner diameter and 1.27 mm wall thickness made of tissue-equivalent plastic TE Shonka 
A-150 (Martino et al. 2010). Under these conditions each TEPC is equivalent to a 
sphere of tissue of 2.75 jum in diameter. 

5. TEPC response to quasi-monoenergetic neutron beams 

It is expected that radiation fields far from the beam axis have a large contribution 
from neutrons produced in fragmentation of beam nuclei. Therefore, before performing 
microdosimetry simulations with nuclear beams we checked the validity of the MCHIT 
model with respect to irradiation of the same TEPC model by neutrons. The 
data obtained with quasi-monoenergetic neutron beams (Nakane & Sakamoto 2001) 
were used for this purpose. Energy spectra of neutrons used in the two sets of 
measurements (Nakane & Sakamoto 2001) were modeled as superpositions of Gaussian 
peaks centered at 40 and 65 MeV, respectively, and broad plateaus. The ratios between 
the flux of neutrons with energies within the peak and the total neutron flux were set to 
1:2.768 and 1:2.807 for 40 and 65 MeV spectra, respectively, following the estimations 
by Nakane et al. (Nakane & Sakamoto 2001). As proved by our simulations, the 
calculated yd(y) distributions are not very sensitive to these ratios. For example, the 
calculations with these ratios set to 1:3.5 provide results similar to ones obtained with 
original ratios (Nakane & Sakamoto 2001). In the following the above- described spectra 
composed of peaks and plateaus are referred to as "40 MeV" and "65 MeV" neutrons. 

MCHIT results for a TEPC located at 5 cm depth inside a PMMA phantom of 
30 x 30 x 30 cm 3 are shown in figure 1 together with the experimental data (Nakane 
& Sakamoto 2001). All microdosimetric distributions for neutron irradiation are 
presented per single energy deposition event in the TEPC. In this experiment the 
phantom was irradiated by "40 MeV" neutrons, and the TEPC was located on 
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the axis of the irradiation field. Neutrons were delivered to the phantom via a 
collimator of 10.9 cm in diameter, which was much larger than the diameter of the 
TEPC. Since neutrons traversed the TEPC randomly at various impact parameters, 
the mean neutron path I inside the TEPC sensitive volume amounted to 2/3 of its 
diameter. As explained in section 2, this corresponds to the standard conditions 
of microdosimetry measurements (International Commission of Radiation Units and 
Measurements (ICRU) 1983). 




10' icr 

y (keV/jim) 



Figure 1. Calculated microdosimetric spectra for a TEPC at 5 cm depth inside a 
PMMA phantom irradiated by quasi-monoenergetic "40 MeV" neutrons. Theoretical 
results obtained with G4QMD and G4BIC models are presented by solid and dashed 
histograms, respectively. Experimental data (Nakane & Sakamoto 2001) are shown by 
points. 



The results presented in figure 1 were obtained by simulating 4 x 10 s primary 
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neutrons for each of the runs where either G4QMD or G4BIC were used to simulate 
neutron-induced nuclear reactions. Microdosimetry spectra calculated for the same 
setup, but for the irradiation with "65 MeV" neutrons are shown in figure 2. As one can 
see from these figures, the general shape of the microdosimetric spectra measured both 
with "40 MeV" and "65 MeV" neutrons is equally well reproduced by MCHIT with 
G4QMD and G4BIC options. At the same time one can note that the results obtained 
with G4BIC are closer to the data for both energies of neutrons. 



10 



n 65 MeV, 5cm 




10 



10' 10* 

y (keV/jim) 



Figure 2. Same as in figure 1, but for quasi-monoenergetic "65 MeV" neutrons. 

Neutrons do not transfer their energy directly to the TEPC gas volume, but rather 
through secondary charged particles produced in neutron interactions with nuclei in 
the gas cavity, TEPC wall or even in surrounding layers of water. Several recoil 
particles can be produced by a single neutron. MCHIT makes it possible to estimate 
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partial contributions from each kind of recoil particles to the total yd(y) curve. Such 
contributions were calculated in the following steps. First, the total dose distribution 
as a function of y was scored in each run together with similar histograms, but scored 
separately for each particle charge. Second, at the end of the run the content of each 
bin of the yd(y) distribution was multiplied by the ratio of the dose from particles of a 
specific charge to the total dose for the same y-bin. The results shown in Figs. 3 and 4 
were obtained in this way. 

Contributions of specific secondary particles to the yd(y) distribution for "65 MeV" 
neutrons are shown in figure 3. These distributions were obtained with G4BIC applied 
to simulate neutron-induced nuclear reactions. This figure demonstrates the origins of a 
peak at y ~ 7 keV//xm and a shoulder extending to higher y values. The peak is due to 
recoil protons, while the shoulder is formed by more heavy recoil nuclei, mostly by alpha- 
particles (y > 30 keV//im) and carbon nuclei (y > 100 keV//im). The contributions 
from boron, lithium and beryllium nuclei are much smaller and attributed mostly to 
y > 200 keV/yum. The microdosimetric distributions calculated for "40 MeV" neutrons 
(not shown) are build in a similar way. However, in this case the maximum is located 
at y ~9 keV//im due to less energetic primary neutrons. Therefore, neutrons transfer 
less energy to recoil protons, which get slower and therefore have higher LET values. 

Contributions of various recoil particles to yd{y) distribution calculated with 
G4QMD for "65 MeV" neutrons are shown in figure 4. As seen from the figure, 
MCHIT predicts more fast recoil protons (with lower LET) compared to experimental 
data. This leads to a shift of the peak of the calculated distribution to lower y. 
Regarding recoil nuclei, results obtained with G4QMD model differ from measurements 
at y ~100 keV//xm where energy deposition events are mostly caused by helium nuclei. 
This can be explained by an underestimation of alpha-particle yields in neutron-induced 
reactions simulated by this model. 

In summary, the MCHIT model with G4BIC applied to simulate neutron- 
induced reactions demonstrates a better agreement with experimental data compared 
to the MCHIT with G4QMD. The maximum of the theoretical yd(y) distributions 
calculated with G4QMD is slightly shifted to lower y values with respect to the 
measurements (Nakane & Sakamoto 2001). However, the overall shape and magnitude 
of the microdosimetric spectra from quasi-monoenergetic neutrons is well reproduced 
also by calculations with G4QMD. The calculated total dose in TEPC increases only 
by 0.3% when G4QMD is replaced by G4BIC. One can note that MCHIT is more 
successful in describing the measured (Nakane & Sakamoto 2001) yd(y) distributions 
as compared with the PHITS code used for modeling (Tsuda et al. 2007) the same 
TEPC measurements (Nakane & Sakamoto 2001). This can be explained by the fact 
that in those PHITS calculations (Tsuda et al. 2007) the production and transport of 
(^-electrons was neglected, while these processes are taken into account in our MCHIT 
simulations. 
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Figure 3. Calculated microdosimetric spectra for a TEPC at 5 cm depth inside a 
PMMA phantom irradiated by quasi-monoenergetic "65 MeV" neutrons. Results of 
MCHIT calculations with G4BIC model are shown. In addition to the total yd(y) 
distribution the contributions of specific particles are shown by various histograms as 
explained on the legend. Experimental data (Nakane & Sakamoto 2001) are shown by 
points. 

6. TEPC response to a therapeutic 12 C beam 

6.1. Number of particles which cross the TEPC 

As described above in section 2, the lineal energy of a particle traversing a spherical 
TEPC can be estimated by dividing the energy e imparted to the detector by the mean 
chord length I = 2 /3d, where d is the TEPC diameter. This estimation is based on the 
assumption that the TEPC is placed in a homogeneous radiation field and traversed by 
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Figure 4. Same as in figure 3, but for MCHIT results with G4QMD model. 



a single particle per event. This corresponds to the simplest standard event topology 
characterized by a single track with I < d. This means that the particle traverses the 
TEPC sensitive volume by a straight-line trajectory, it is not stopped inside the volume, 
and secondary particles are not produced neither in the volume nor in the vicinity of it. 

In the experiment (Martino et al. 2010) the TEPC was irradiated either directly by a 
focused 12 C beam, or hit only by secondary particles at TEPC locations outside the beam 
spot. In this section we investigate the event topology relevant to the both cases. Indeed, 
on the contrary to the measurements with neutron beams described in section 5, the 
size of the 12 C beam was essentially smaller (3 mm FWHM) compared to the diameter 
of the TEPC sensitive volume (d = 12.7 mm). Therefore, for a TEPC placed on the 
beam axis at the plateau region of the depth-dose distribution and, possibly, close to the 
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Bragg peak, the mean chord length will be larger than 2 /3d, as small impact parameters 
dominate. Moreover, several secondary particles can impact the TEPC in a single event. 
In this section we investigate the influence of these conditions on the determination of 
y f and y d . 

Since a TEPC detector is only sensitive to the total energy e deposited to its volume 
given by the number of ionized atom-electron pairs produced in the gas chamber, the 
number of tracks in each event and their length (event topology) can not be identified 
in experiment. However, one can retrieve this information from MCHIT simulations by 
scoring the number of tracks and their total length inside the TEPCs. The probability 
distributions for the number of tracks f(N) inside the TEPCs are shown in figure 5. 
The distributions are calculated for the events with at least one particle traversing the 
TEPC, and they are given at nine positions in the water phantom listed in table 2. These 
distributions are calculated per single beam particle and represent the probability for 
a TEPC placed at a certain position in the phantom to be traversed by N particles 
of any kind. This number is calculated by means of the Geant4 primitive scorer 
G4PSCellCurrent and it amounts to the number of particle tracks N which pass through 
the TEPC, but do not start or stop inside the TEPC. In particular, this means that 
numerous low-energy ^-electrons which are produced by energetic nuclei and stopped 
inside the TE gas are not counted by this scorer. Following this definition of the scorer, 
N — 1 corresponds to the standard event topology defined above despite of multiple 
(^-electrons stopped inside the TEPC. The tracks of the particles which do not deposit 
energy to the TEPC are also included. 

The distributions presented in figure 5 were calculated with and without production 
of (^-electrons. In the latter case only secondary nucleons and nuclear fragments were 
produced and counted while they traversed TEPCs. As one can see from the distribution 
calculated at "0 cm, plateau", the TEPC at this location is traversed by beam particles 
with ~ 0.85 probability. This reflects the attenuation of the primary beam in the first 
5 cm of water. On the contrary, the probability to hit the TEPC at "10 cm, tail" 
position far from the beam axis is only ~ 0.0053 per beam particle. As deduced from 
the comparison of the distributions calculated with and without generating (^-electrons, 
all TEPCs placed at the beam axis are usually traversed by several electrons in addition 
to nucleons and nuclei. 

The topology of each event associated with the traversing of TEPC by one or several 
particles, can be characterized by the average number of tracks N calculated per event. 
These average numbers are presented in each panel of figure 5 for both calculational 
options - with and without (^-electrons. When the production of 5-electrons is neglected, 
the most probable events are characterized by a single particle track, with N — 1.5 and 
1.4 under the direct impact of the beam at the plateau and peak, respectively. However, 
when the production of (^-electrons is considered, N = 3.6 at the plateau and N = 2.7 at 
the Bragg peak position. Although these electrons deliver less energy to the TEPC when 
compared to beam nuclei, their contribution can affect the results of simulations. As 
demonstrated in figure 5, the A-distributions and N calculated for other seven positions 
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Figure 5. Probability distributions for the number of tracks in TEPCs at different 
locations in the water phantom irradiated by 300 A MeV 12 C beam. The condition 
to have at least one track in the corresponding TEPC is applied. The distributions 
with and without simulation of (5-electrons are presented by solid- and dashed-linc 
histograms, respectively, and calculated per beam particle. The average number of 
tracks N per traversing event are presented for each TEPC location, and the values 
obtained without 5-electrons are given in parentheses. 



are less affected by neglecting 5-electrons. At these positions the TEPCs are not directly 
impacted by the 12 C beam. 

The probability distributions f(l) for the total track length I calculated with 
G^PSPassageTrackLength for particles traversing the TEPC at different locations in 
the water phantom are shown in Figs. 6 and 7. The values of l/d calculated per 
event of TEPC traversing are also given for each TEPC position. We first consider 
the distributions obtained without production of 5-electrons, figure 6. In this case all 
tracks simulated in each TEPC are only due to nucleons or nuclei. 

As expected, the irradiation of TEPCs by a thin 12 C beam is reflected in the 
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Figure 6. Probability distributions for the total track length for particles traversing 
the TEPCs at different locations in the water phantom irradiated by 300A MeV 
12 C beam. The condition to have at least one track in the TEPC is applied. The 
calculations were performed without production of (5-electrons. The average values of 
l/d calculated per event are listed for each TEPC position. 



/(^-distributions calculated for "0 cm, plateau" and "0 cm, peak". The values of 
l/d = 1.284 and 1.129 in these cases are about twice as large compared to l/d = 0.667 
for homogeneous radiation field. The linear dependence of /(/) versus / (for / < d) 
indicates a random particle incidence for the TEPC placed on the beam axis, but beyond 
the Bragg peak ("0 cm, tail"). Nevertheless, due to multiple tracks in a single event 
l/d is still large, l/d ~ 0.8. Some events with two and three nucleons or nuclei in the 
TEPC are characterized by I > d at "0 cm" positions, as already demonstrated earlier 
in figure 5. The distributions at the radius of 2 cm still exhibit the presence of one and 
two tracks of nucleons or nuclei in the TEPC. The mean chord length is close to 2 /3d 
only for the TEPCs at the radius 10 cm. There the probability of an event with two 
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particles in the TEPC is negligible. 

The f(l) distributions calculated with production of (^-electrons are shown in 
figure 7. As follows from the comparison of Figs. 6 and 7, secondary electrons noticeably 
change the f(l) distributions only in the cases when the TEPC is directly irradiated by 
the 12 C beam, at "0 cm, plateau" and "0 cm, peak" positions. This comparison also 
demonstrates that the events with / > 2d and, respectively, with N > 2 (figure 5) are 
mainly associated with the production of 5-electrons rather than nucleons or nuclear 
fragments. 




Figure 7. Same as in figure 6 but simulated with production of (5-electrons. 

In summary, the analysis of f(N) and f(l) distributions demonstrates different 
topologies of events in the TEPCs placed at the beam axis and far from the beam. In 
the latter case the events are mostly characterized by the simplest standard topology 
associated with a single track of nucleon or nuclear fragment traversing the TEPC 
sensitive volume. On the contrary, the events in the TEPCs placed on the beam axis at 
the entrance to the phantom ("0 cm, plateau") and at the Bragg peak ("0 cm, peak") 
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are more complicated. They are characterized by one or two tracks of beam nuclei, 
secondary nucleons and/or nuclear fragments traversing the TEPC sensitive volume. In 
addition, tracks of two, three or more energetic electrons are frequently present in the 
TEPCs placed on the beam axis. 

6.2. Lineal energy spectra inside the water phantom 

Following the consideration given above to particle track patterns in TEPCs, we consider 
now the corresponding calculated microdosimetric yd(y) distributions presented in 
Figs. 8 and 9. As discussed in section 6.1, for the TEPCs placed on the beam axis, 
I > 2 /3d since a certain number of events is characterized by more than one track 
traversing the TEPC. However, in GSI measurements (Martino et al. 2010) I — 2/3d was 
assumed following the standard microdosimetry technique (International Commission of 
Radiation Units and Measurements (ICRU) 1983). The same normalization was taken 
in our calculations of y for the sake of a direct comparison with GSI data. 

The MCHIT calculations were performed with and without modeling secondary 
electrons, and also with two different nuclear fragmentation models. The sensitivity 
of calculated distributions to computational parameters can be studied by comparing 
these distributions with each other and with experimental data. 

The microdosimetric yd(y) distributions measured at GSI (Martino et al. 2010) 
inside the water phantom irradiated by a 300A MeV 12 C beam are shown in figure 8 
together with our simulation results J. The distributions are given per beam particle 
for the same TEPC positions inside and outside the beam spot as listed in table 2 and 
described in section 6.1. In order to reach a better agreement between measured and 
calculated positions of the peak seen in yd(y) distributions at "0 cm, peak", in our 
simulations the TEPCs were shifted by 2 mm deeper in water. One can explain this 
shift by uncertainties in the measurements of either TEPC positions or beam energy, as 
well as by some inaccuracy of the electromagnetic physics models of Geant4. 

The importance of 5-electrons can be estimated by considering yd(y) distributions 
obtained with and without simulating electron production, but with the same nuclear 
fragmentation model (G4QMD), as shown in figure 8. Here, the distributions calculated 
with these two options noticeably differ at the "0 cm, plateau" TEPC position close to 
the beam entrance to the phantom. In particular, the number of energy deposition 
events with 20 < y < 200 keV//zm is much higher compared to the calculation with 
^-electrons. At this point fast beam nuclei produce energetic electrons which may escape 
from the TEPC sensitive volume and thus reduce the energy deposited to this volume. 
However, the yd(y) distribution calculated at the Bragg peak (at "0 cm, peak") is not 
affected by neglecting ^-electrons. Much more secondary electrons are produced close 
to the Bragg peak, but they are generally less energetic compared to those at "0 cm, 

| On the contrary to the standard definitions of yd(y) spectra given in section 2, the data (Martino 
et al. 2010) were not divided by yj. Therefore, this factor was also omitted in our simulation results 
presented in Figs. 8, 9 and 12. 
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plateau" and do not escape from the TEPC. This means that the energy deposited to the 
TEPC at "0 cm, peak" can be simply calculated from the stopping power of beam nuclei 
even without modeling (^-electrons. As seen from figure 8, the production of (^-electrons 
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Figure 8. Microdosimetric yd(y) distributions in water phantom irradiated by 
300 A MeV 12 C nuclei. Distributions calculated by MCHIT per beam particle with 
and without simulating ^-electrons are presented by solid- and dashed-line histograms, 
respectively. Points represent experimental data (Martino et al. 2010). 

also changes the distributions in the range of 10 < y < 100 keV/^m for the TEPCs 
located at 2 cm from the beam axis. Similar, but smaller changes in yd(y) distributions 
associated with electron production are also found at 10 cm radii. At all these positions 
energy is deposited to the TEPC mostly by protons, and their contribution is reduced 
if secondary electrons are produced and escape the TEPC volume. 

MCHIT results obtained with G4BIC and G4QMD models used to simulate nuclear 
fragmentation are presented in figure 9 together with experimental data (Martino 
et al. 2010). One has to keep in mind the range of magnitudes of the distributions 
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presented in nine panels of Figs. 8 and 9. They are normalized per single beam 
particle and have a span of almost six orders of magnitude. This is because of the 
fact that TEPC hits at the position "10 cm, plateau" are approximately 10 4 times less 
frequent than at the Bragg peak ("0 cm, peak"), as described below in section 6.3. The 
TEPCs located at 10 cm distance from the beam axis are hit exclusively by secondary 
particles (mostly nucleons) produced in nuclear fragmentation reactions. Therefore, the 
microdosimetry data (Martino et al. 2010) provide another possibility to validate the 
nuclear fragmentation models of Geant4. 

The experimental data (Martino et al. 2010) were reported only above 0.3- 
0.5 keV//xm, while the ^(^-distributions were calculated without restricting y values. 
In this way MCHIT predicts the microdosimetric spectra also in the domains where 
data are not available. 

Detailed consideration of the contributions of different charged particles to the 
calculated microdosimetric spectra will be given elsewhere (Burigo et al. 2012). In this 
work only general features of yd(y) distributions are discussed, and the attention is 
focused on their general shape due to primary nuclei or secondary particles. 

In the measurements (Martino et al. 2010) a sharp peak in the spectrum is observed 
at the TEPC position "0 cm, plateau", figure 9. It is located at y ~16 keV//zm, and its 
position is well reproduced by the MCHIT simulations with both G4BIC and G4QMD. 
However, as seen from the corresponding panel of figure 9, MCHIT predicts a slightly 
sharper drop of yd(y) at the right slope of the peak. This peak is due to energetic 
beam nuclei that traverse the TEPC at the entrance to the phantom. As discussed in 
section 6.1, because of a relatively small beam diameter (3 mm FWHM) compared to 
the TEPC size, carbon nuclei propagate close to the diameter of the TEPC. As seen in 
the panel "0 cm, peak", the peak in the yd(y) distribution corresponding to 12 C nuclei 
becomes higher and broader as the TEPC is moved to the stopping point of beam 
nuclei. Compared to the "0 cm, plateau" panel, here the peak is shifted to a much 
higher value of y ~131 keV//im, and it has a prominent satellite peak with a maximum 
at y ~25 keV//xm. The positions of both peaks are well reproduced by simulations. The 
observed shift of the main peak to larger y is due to increasing the LET of beam nuclei 
as they are slowing down with their penetration in water. 

A common feature of the "0 cm, plateau" and "0 cm, peak" distributions consists in 
the presence of additional peaks at lower y values compared to the peaks of beam nuclei. 
Such satellite peaks are present due to relatively heavy projectile fragments (boron, 
beryllium and lithium nuclei) propagating with velocities close to the velocity of beam 
nuclei. Since their Z 2 are smaller, their ionization energy loss is reduced accordingly. 
This leads to reduced LET of secondary fragments compared to beam nuclei. While 
the satellite peak is accurately reproduced by theory at the entrance to the phantom 
at "0 cm, plateau", a similar peak at "0 cm, peak" TEPC position is underestimated. 
This indicates, that the yields of secondary fragments are underestimated by the nuclear 
fragmentation models used in MCHIT. 

The spectrum at the beam axis beyond the Bragg peak ( "0 cm, tail" ) and all the 
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Figure 9. Microdosimetry spectra in water phantom irradiated by 12 C 300 A MeV 
nuclei calculated by MCHIT with G4QMD and G4BIC options for nuclear 
fragmentation, shown as solid- and dashed-line histograms, respectively. Points 
represent experimental data (Martino et al. 2010). 



six spectra at 2 and 10 cm radii do not demonstrate any sharp peaks. This is because 
of contributions from various secondary light particles characterized by a broad range 
of kinetic energies and charges. It is expected that recoil charged particles produced 
in neutron-induced reactions also contribute to these spectra. This contribution is 
estimated below in section 7. While the general trends of the distributions measured at 
these seven points are reproduced by MCHIT, the absolute values are underestimated 
by a factor of two or slightly less. In view of a 10 4 difference in magnitude between 
"0 cm, tail" and "10 cm, plateau" distributions, such level of agreement between MCHIT 
and data can be accepted. Nevertheless, this discrepancy suggests that there is still a 
room for improvements of the considered nucleus-nucleus collision models of Geant4, 
G4BIC and G4QMD, or de-excitation models. 
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6.3. Hit probability, yf, and dose inside the water phantom 

We consider now integral microdosimetric quantities yf, yd defined in section 2 as well 
as the total dose at various positions inside the water phantom irradiated by 300A MeV 
12 C nuclei. The total dose calculated at various distances from the beam axis is shown in 
figure 10 as a function of the depth in water together with measured dose values (Martino 
et al. 2010). 
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Figure 10. Calculated and measured total doses at various depths in the water 
phantom irradiated by 300^4 MeV 12 C nuclei. The doses calculated at 0, 20 and 
100 mm radii by scoring energy imparted to the TEPC gas cavity are shown by open 
squares. The doses calculated at the same positions from microdosimetric distributions 
according to Eq. (4) are shown by open triangles. Lines of various styles explained on 
the legend present total dose calculated in water by direct scoring of energy deposition 
at 0, 10, 20, 50 and 100 mm radii. The measured doses (Martino et al. 2010) are shown 
by open circles connected by dark solid lines to guide the eye. 



As explained by the legend of figure 10, calculated doses were obtained by three 
different methods. In the first method labeled as "TEPC" the dose delivered to the 
tissue-equivalent gas was directly scored in simulations. For this purpose the TEPC gas 
cavity was defined as a sensitive volume, and a primitive scorer GJ^PSEnergyDeposit was 
used. In the second method labeled as "micro-calc" the dose value was calculated from 
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microdosimetric distributions following Eq. (4). Finally, in the third method calculations 
were performed in water without placing TEPC. In the latter method a set of concentric 
rings was defined inside the phantom with 0.1 mm steps in the depth and radius. At the 
end of each run the energy imparted to each ring was calculated and then divided by the 
mass of the ring and by the number of beam particles. Due to a finite size of the TEPC, 
particles traverse it at various distances from the beam axis. Therefore, in calculating 
the dose without TEPC one has to take an average dose value for a characteristic volume 
of a similar size. For example, the average dose for 0-6.35 mm radii was calculated to 
obtain the dose on the beam axis. The radial ring thickness for dose calculations was 
set to 2/3 of the TEPC diameter at other distances from the beam axis. 

The calculated dose values obtained by "TEPC" and "micro-calc" methods agree 
well with each other and with the data at most of the TEPC positions, see figure 10. 
At the same time, the dose at "0 cm, plateau" is underestimated, as one could expect 
from microdosimetric spectra calculated at this TEPC position, see section 6.2. The 
dose values calculated in the water phantom without TEPC agree well with the TEPC- 
based measurements and TEPC-based simulations except the points in the vicinity of 
the Bragg peak. As expected, these points are characterized by high dose gradients 
which are predicted by calculations without TEPC, e.g. at 10 mm radius. These dose 
gradients are smoothed by taking average values in water over the above- described rings. 
However, the dose in TEPC-based measurements still differs from the dose calculated 
without TEPC in the Bragg peak region at 10 mm radius. 

Figure 11 gives further insight into the influence of a finite size of TEPC on 
microdosimetry measurements with therapeutic pencil-like beams. The doses from 
TEPC-based simulations, calculations without TEPC and measured doses agree well 
with each other at 10 and 20 mm distance from the beam axis. The same is true for the 
measurements performed at the beam axis in the tail. In all these cases the spatial dose 
distribution is characterized by relatively small gradients inside TEPC volumes which 
suggests good agreement between calculations with and without TEPC. While a good 
agreement is observed between data and TEPC-based simulations, direct calculations 
of the dose differ from both of them on the beam axis at the plateau and peak. As 
seen, the local dose at the TEPC center is generally 2-3 times higher than the TEPC- 
measured and TEPC-simulated dose under the focused irradiation of this TEPC by 
3 mm FWHM pencil-like beam. Our results suggest that a good agreement between 
measured and calculated doses can be obtained only by direct modeling of the TEPC 
geometry, and not by scoring the dose in simulations without TEPC, as done by other 
authors (Hultqvist et al. 2010, Taleei et al. 2011). 

The calculated frequency-mean yf and dose- mean lineal energies are given in 
tables 3 and 4. The corresponding values calculated from experimental microdosimetric 
distributions are also given in these tables for comparison. The values presented in 
table 3 are given per beam particle, and through Eq. (4) they are connected to the 
total dose at various locations inside the phantom. In addition, the probability Ptepc 
of energy deposition at a given TEPC position is also listed in table 3. From these 
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Figure 11. Calculated and measured total dose in the water phantom irradiated by 
12 C nuclei as a function of the distance from the beam axis. The doses calculated on 
the beam axis and 20 mm away by scoring energy imparted to the TEPC gas cavity 
are shown by open squares. Lines of various styles explained on the legend present 
total dose calculated in water by direct scoring of energy deposition at the plateau 
(44.04 mm depth), Bragg peak (171.04 mm) and tail (269.04 mm). The measured 
doses (Martino et al. 2010) are shown by open circles. Lines with arrows demonstrate 
the radial ranges covered by the TEPC placed at 0, 10 and 20 mm from the beam axis. 

probability values one can conclude, for example, that on average there is one deposition 
event per ~ 2000 beam particles for the TEPCs placed 10 cm away from the beam axis, 
both at the peak and tail. The MCHIT results for yj obtained with the standard 
electromagnetic model and with the Penelope model differ by ~ 5 — 10% at all TEPC 
positions. The measured y~f values on the beam axis at the plateau and peak are well 
described by calculations. However, the theoretical values are lower by a factor of 1.5 — 2 
compared to the measured ones at the other seven TEPC positions, where TEPCs are 
mostly hit by secondary protons and neutrons. 

The calculated and measured values are presented in table. 4. They are given 
per energy deposition event according to Eq. (2) (note 1/jjf factor in this equation), 
and related to the radiation quality at the corresponding TEPC locations in the 
phantom. With the exception of "10 cm, tail" position, y d calculated with the standard 
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Table 3. Probabilities to deposit energy in the TEPC at various positions and the 
corresponding yf, all calculated per beam particle. The calculations were performed 
with Standard Electromagnetic Physics (Option 3) physics list. The results obtained 
with G4EmPenelope model for transporting electrons, positrons and gammas are 
given in the third column. The values of yf calculated from experimental (Martino 
et al. 2010) distributions yd(y) are given for comparison. 



TEPC's 
position 


-Ptepc 

MCHIT 
Standard opt3 


Vf 

(keV/yum) 
MCHIT 
Standard opt3 


Vf 

(keV/Mm) 

MCHIT 

Penelope 


Vf 

(keV/yum) 

GSI experiment 


cm, plateau 


9.88- 10- 1 


14.6 


14.9 


16.9 


cm, peak 


6.53- 10- 1 


63.7 


65.3 


71.2 


cm, tail 


4.15- 10~ 2 


1.44- lO" 1 


1.59- 10- 1 


3.19- 10- 1 


2 cm, plateau 


5.67- 10- 3 


9.29- 10~ 3 


9.83- 10~ 3 


1.93- 10~ 2 


2 cm, peak 


1.72- 10~ 2 


3.09-10- 2 


3.39- 10~ 2 


5.83- 10~ 2 


2 cm, tail 


1.58- 10~ 2 


4.12 • 10~ 2 


4.42 • 10~ 2 


7.42 • 10~ 2 


10 cm, plateau 


3.03 • 10~ 5 


7.07- 10" 5 


6.96 • 10~ 5 


1.34- 10~ 4 


10 cm, peak 


2.66- 10~ 4 


4.45 • 10~ 4 


4.80 • 10~ 4 


7.27- 10~ 4 


10 cm, tail 


2.80 • 10~ 4 


4.63 • 10~ 4 


5.02 • 10~ 4 


7.49 • 10~ 4 



electromagnetic and Penelope models agree well with each other. Since yd values were 
obtained by integrating y 2 f(y)/yf distributions instead of yf(y) for yf, the results are 
less sensitive to low-energy thresholds set in measurements. This explains the improved 
agreement of the calculated ya with data as compared to £//. As seen from the table, 
calculated with the standard electromagnetic and Penelope models at the beam axis and 
at "2 cm, peak" and "10 cm, peak" locations agree well with data. However, differences 
between MCHIT results and data are seen at "2 cm, plateau", "10 cm, plateau" and 
"10 cm, tail". As only protons and neutrons contribute at these TEPC positions far 
from the beam, this discrepancy points at the necessity to improve the description of 
production of these secondary particles by MCHIT. 

6.4- Relations between yf and LET for various beam profiles and surrounding media 

As discussed in section 6.1, the TEPC response to focused and homogeneous irradiation 
is different. It can be also influenced by the media which surrounds the TEPC, as 
particles produced in this media, e.g. in fragmentation reactions, also hit the detector. 
It is expected (International Commission of Radiation Units and Measurements 
(ICRU) 1983) that the relation: 

Vf = L, (5) 

should hold for a spherical TEPC randomly traversed by particles with a constant 
LET L, in the case when the production of (^-electrons and other secondary particles 
is neglected. The relations between measured, simulated y~f and LET under various 
conditions can be assessed by considering Tab. 5. In addition to y/ calculated with 



Microdosimetry of radiation fields from 12 C beams 



30 



Table 4. Calculated and measured (Martino et al. 2010) yd per energy deposition 
event in each TEPC. The calculations were performed with Standard Electromagnetic 
Physics (Option 3) physics list. The results obtained with GJ^EmPenelope model for 
transporting electrons, positrons and gammas are given in the second column. The 
experimental values were obtained by integrating yd(y) distributions and dividing the 
results by the experimental yf values listed in table 3. 



TEPC's position 


yd 


yd 


yd 




(kcV/Mm) 


(keV//xm) 


(keV/^tm) 




MCHIT 


MCHIT 


GSI experiment 




Standard opt 3 


Penelope 




cm, plateau 


17.0 


17.3 


17.2 


cm, peak 


188. 


192. 


162. 


cm, tail 


14.9 


16.0 


13.6 


2 cm, plateau 


13.8 


9.25 


7.05 


2 cm, peak 


9.39 


11.2 


8.64 


2 cm, tail 


12.4 


13.1 


9.34 


10 cm, plateau 


17.9 


17.6 


12.9 


10 cm, peak 


12.5 


11.1 


10.1 


10 cm, tail 


12.3 


20.4 


8.83 



default physics settings of MCHIT and labeled with (a); the results obtained with the 
Penelope model used for transporting electrons (b); without simulating the production 
and transport of (^-electrons (c); and without simulating nuclear reactions (d) are 
presented. 

In Tab. 5 the results for the TEPC placed at "0 cm, plateau" at the depth of 
52.1 mm in water under the impact of 3 mm FWHM beam are presented, cases (la), (lb), 
(lc) and (Id). At this TEPC position the energy of 12 C at the entrance to the TEPC 
volume is estimated as ~218A MeV, and it is close to the beam energy (220 A MeV) 
used by Taddei et al. for homogeneous irradiation of a similar TEPC§, but surrounded 
by air (Taddei et al. 2008). This makes possible the comparison of various irradiation 
conditions presented in table 5. 

The y f calculated with MCHIT for TEPC in water 14.6 keV//im (la) is relatively 
close to LET, but remains smaller compared with the experimental result of GSI 
(16.9 keV//xm). In the case (lb) the Penelope model with the production and transport 
of ^-electrons extended to lower electron energy provides a larger yf, of 14.9 keV/ jum, due 
to the enhancement of energy deposition inside the gas cavity by low-energy (^-electrons 
produced in the plastic shell. In the case (lc) the production of (^-electrons is neglected in 
the calculation, and y~f increases to 16.1 keV//im and exceeds the corresponding LET. 
This is because of additional energy deposition in the TEPC sensitive volume which 

§ The thickness of the TEPC wall was twice as large (2.54 mm) compared to the TEPC used at 
GSI, and their TEPC operated at different tissue-equivalent gas pressure. Furthermore, those energy 
deposition events in the TEPC due to outcoming nuclei different from the incoming beam nuclei were 
ruled out in the experimental method. All these details were taken into account in our simulations. 
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otherwise would be taken away by secondary electrons propagating beyond the TEPC. 
The case (Id) presents calculations performed with taking into account secondary 
electrons, but neglecting beam fragmentation. This means that instead of the mixture 
of 12 C beam nuclei and secondary fragments, which hit the TEPC in the cases (la), 
(lb) and (lc), the detector is traversed only by 12 C nuclei and due to their higher Z 2 
factor, i/f increases to 16.5 keV//xm, well above LET. This also explains the value of 

Table 5. Calculated and measured frequency-mean lineal energy jjf per energy 
deposition event in the TEPC. Results are given for a Gaussian 3 mm FWHM beam 
profile (1,2) and for homogeneous irradiation of the TEPC (3,4) surrounded by water 
(1) or air (2,3,4). Calculations with various physics settings are marked by letters: (a) - 
default MCHIT physics, (b) - the Penelope model for electrons, (c) - without simulation 
of 5-electrons, (d) - without simulation of nuclear fragmentation. Experimental 
data (Martino et al. 2010), (1) and (Taddei ct al. 2008), (3,4) are given for comparison 
with MCHIT results. LET values calculated (Toftegaard & Bassler 2012, Luhr 
et al. 2012) according to ICRU Report 73 (International Commission of Radiation 
Units and Measurements (ICRU) 2005) are also given. 



case 


beam 


beam 


media 


Vf 


Vf 


Vf 


LET 




energy 


profile 






Calculation 


Experi- 


ICRU 














ment 


73 










MCHIT 


Taddei 


Martino 














et al. 


et al. 
















Taddei 






(MeV/ 










et al. 






nucleon) 






(keV//xm) 


(keV//xm) 


(keV//xm) 


(keV//xm) 


(la) 


218 


Gaussian 


water 


14.6 




16.9 


15.23 


(lb) 


218 


Gaussian 


water 


14.9 




16.9 


15.23 


(lc) 


218 


Gaussian 


water 


16.1 




16.9 


15.23 


(Id) 


218 


Gaussian 


water 


16.5 




16.9 


15.23 


(2a) 


220 


Pencil-like 


air 


20.4 






15.14 


(3a, d) 


220 


flat 


air 


13.8 


14.5 


13.4 


15.14 


(3b,d) 


220 


flat 


air 


14.1 


14.5 


13.4 


15.14 


(4a,d) 


390 


flat 


air 


9.3 


9.4 


9.93 


11.02 


(4b,d) 


390 


flat 


air 


9.99 


9.4 


9.93 


11.02 



yf calculated for the conditions of the second experiment, case (2). There the beam 
fragmentation outside the TEPC is negligible as it is surrounded by air. The TEPC is 
traversed almost exclusively by 12 C nuclei along the TEPC diameter. This explains the 
largest energy deposition, and hence largest y~f, obtained for this case. 
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In order to demonstrate the dependence of y~f on the beam shape, the case of a 
flat beam (3,4) was also considered. In this case y~f = 13.8 keV//im which is noticeably 
smaller than y~f of the case (2), 20.4 keV//xm, calculated for the pencil-like beam of the 
same energy. The ratio between y~f for the random impact on the TEPC (3) and for the 
central beam incidence (2) is remarkably close to the expected value of 2/3. One can 
note that the MCHIT result for the flat beam is closer to the experimental value than 
the value calculated by Taddei et al. with another version (7.1) of the Geant4 toolkit. 
We attribute this improvement to the changes in electromagnetic models done since the 
release of the version 7.1 and to differences in calculational parameters. 

Finally, the cases (4a,d) and (4b, d) represent calculated and measured yf for the 
TEPC in air under a random impact of a 390 A MeV beam (Taddei et al. 2008). The 
y~f value calculated with the standard electromagnetic models is lower compared to the 
measured £// (Taddei et al. 2008). All the calculated values for the case (4) remain 
below LET. The yf calculated by MCHIT with the Penelope model (4b, d) agrees better 
with the experimental value at 390A MeV beam energy. 

Depending on the beam profile the measured y~f are either ~ 20% higher than LET 
for focused irradiation, case (1), or ~ 10% lower than LET for homogeneous irradiation, 
case (3). As follows from this analysis, the exact modeling of TEPC irradiation 
conditions (beam profile and surrounding media) is important for reproducing data. 
As demonstrated by calculations, the simulation results are rather insensitive to the 
choice of electromagnetic models. This shows that the production and transport of 
^-electrons between 100 and 990 eV may not play an important role for the energy 
deposition inside the gas cavity. 

7. Contributions of secondary neutrons to microdosimetric quantities 

The assignment of energy-dependent weighting factors is critical for calculating 
equivalent doses for neutrons, e.g. in radiation protection. As discussed in several 
publications, the calculations of dose from secondary neutrons in proton (Jiang 
et al. 2005, Jarlskog & Paganetti 2008, Xu et al. 2008) and heavy-ion (Newhauser 
& Durante 2011) therapy are prone to various uncertainties. First, there are common 
difficulties for all kinds of radiation in estimating the risks from low doses (Kellerer 2000). 
Second, depending on the kind of tissue under irradiation in vivo and neutron energy, 
the relative biological efficiency (RBE) varies from 2 to 50 (Grahn et al. 1992) or even 
from 7 to 70 (Dennis 1987). Third, the fraction of neutrons in complex radiation fields 
surrounding therapeutic beams should be also properly evaluated, as discussed in this 
section. 

The calculations presented above in section 5 demonstrate that fast neutrons 
deposit energy to the tissue-like media mostly by recoil protons. Neutrons produced 
in fragmentation reactions during particle therapy can travel long distances before they 
interact, thus depositing energy all over the patient's body. This issue was specially 
investigated with the MCHIT (Pshenichnov et al. 2005) with the conclusion that the 



Microdosimetry of radiation fields from 12 C beams 



33 



neutron doses in carbon and proton therapy are at the same order, about 1% of the 
total dose. As demonstrated by recent measurements (Yonai et al. 2010), the dose 
from secondary neutrons in passive carbon-ion therapy is indeed comparable to the 
neutron dose in proton therapy. However, it is not yet clear whether neutrons produced 
during proton or carbon treatments essentially elevate the risk of secondary cancer (Jiang 
et al. 2005, Jarlskog & Paganetti 2008, Xu et al. 2008, Newhauser & Durante 2011). 
The contributions of secondary neutrons to the microdosimetric spectrum can be 




y (keV/(im) y (keV/p,m) y (keV/p.m) 

Figure 12. Microdosimetry spectra in water phantom irradiated by 12 C 300^4 MeV 
nuclei calculated with MCHIT. Simulations are performed with G4QMD model. Points 
represent experimental data (Martino et al. 2010). 

measured applying a veto counter technique which classify events due to charged or 
neutral particles. This technique was used elsewhere (Endo, Tanaka, Takada, Onizuka, 
Miyahara, Sato, Ishikawa, Maeda, Hayabuchi, Shizuma & Hoshi 2007, Wissmann 
et al. 2010) for 290^4 MeV carbon beam in acrylic phantom and 200A MeV carbon 
beam in water phantom, respectively. Endo et al. found that the neutron contribution 
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to the deposited dose in the tail was 18% for the beam axis and 51% at 10 cm radius. 
For larger distances from the beam axis the neutron dose is predominant, increasing 
from 74% in forward direction to 89% in backward direction at 15 cm radius. 



Table 6. Dose inside the TEPC due to secondary neutrons. The column "TEPC" 
stands for the dose obtained by scoring the energy deposited inside the TEPC. The 
column "micro-calc" stands for the dose calculated from the mean lineal energy. 



TEPC's 
position 
notation 



Total Dose (nGy/ion) Neutrons' Dose (nGy/ion) 

GSI experiment TEPC micro-calc TEPC % micro-calc % 



cm, 
plateau 
cm, 
peak 
cm, 
tail 
2 cm, 
plateau 
2 cm, 
peak 
2 cm, 
tail 
10 cm, 
plateau 
10 cm, 
peak 
10 cm, 
tail 



30.6 
93.4 

2.74 • 10" 1 
1.68 ■ 10- 2 
5.04 • 10~ 2 
6.42 • 10~ 2 
1.20 • 10~ 4 
6.25 • 10" 4 
6.44 • 10" 4 



18.4 



81.3 



18.4 



80.6 



1.88 -lO" 1 1.88 • 10" 1 



1.43-10" 2 1.41-10" 



4.50 • 10~ 2 4.52 • 10" 2 



5.63 • 10~ 2 5.56 • 10" 2 



9.17 -10- 5 9.83 -10" 5 



6.33 • 10~ 4 6.49 • 10" 4 



7.16 • 10" 4 6.81 • 10" 4 



< 5.0 • 10- 3 < 0.03 < 5.0 • 10~ 3 < 0.03 



< 5.0 • 10~ 2 < 0.06 < 5.0 • 10~ 2 < 0.06 



2.9 • 10~ 3 1.5 1.1 • 10~ 3 0.6 



9.2 • 10" 4 6.4 7.8 • 10" 4 5.5 



1.5 • 10" 3 3.3 2.0 • 1(T 3 4.4 



2.0 -10~ 3 3.6 1.6 • 10 



3.7 -10- 5 40 



1.6 • 10~ 4 23 



1.8 • 10~ 4 25 



2.9 



4.3 • 10" 5 44 



1.5 • 10" 4 23 



1.8 • 10~ 4 26 



In GSI measurements (Martino et al. 2010) contributions to microdosimetric 
quantities from secondary neutrons were not identified. Monte Carlo simulations with 
MCHIT make it possible to estimate the energy deposition due to interactions of 
neutrons. This is performed in two steps. Firstly, the energy deposition and lineal 
energy spectra are calculated taking into account all relevant physics processes, including 
production of secondary neutrons and their interactions with the phantom and TEPC. 
Secondly, calculations were performed for the same set-up including neutron production, 
but neglecting secondary interactions of neutrons. Thus, the difference between these 
two simulations provides the neutron contributions to the dose and yd(y) distributions. 

Two different methods were used to calculate the dose in TEPC due to secondary 
neutrons. In the first method the dose inside the TEPC, which is represented 
by a corresponding sensitive volume, is scored directly. Then neutron contribution 
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is estimated as a difference between two calculations with and without secondary 
interactions of neutrons. In the second method the dose is estimated from y~f calculated 
from microdosimetric spectra of neutrons. The absolute dose values obtained by these 
two methods and marked as "TEPC" and "micro-calc" , respectively, are given in table 6. 

The contribution from secondary neutrons to the microdosimetric spectra is shown 
in figure 12. As seen, neutron contribution increases with the distance from the beam 
axis. Both calculational procedures show good agreement in estimating neutron dose. 
The neutron contribution at 10 cm radius increases from ~ 25% in the forward direction 
to ~ 40% in the backward direction. These results can not be directly compared to those 
of Endo et al. (Endo, Tanaka, Takada, Onizuka, Miyahara, Sato, Ishikawa, Maeda, 
Hayabuchi, Shizuma & Hoshi 2007) as they correspond to different TEPC position and 
beam profiles. However, the observed trend of increasing neutron contribution with 
respect to the distance from the primary beam is in agreement with their experimental 
findings. One should bear in mind, however, that the absolute doses are very small in 
these regions. 

8. Conclusions 

Since a TEPC placed in a water phantom irradiated by nuclear beam inevitably changes 
the amount of material which is traversed by particles, the radiation field is distorted in 
the presence of the detector. The TEPC geometry was thoroughly implemented in the 
Monte Carlo simulations performed in this work with the MCHIT model. This allowed 
to obtain a good agreement between calculated and measured yd(y) distributions, £//, 
y d and the total dose for different TEPC positions inside the water phantom. 

Our calculations of the microdosimetric variables and the total dose in the water 
phantom irradiated by the pencil-like beam of 300 A MeV 12 C nuclei lead us to the 
following conclusions. 

• Contributions of primary beam nuclei and secondary fragments can be distinguished 
in the calculated and measured yd(y) microdosimetric spectra on the beam axis. 

• As follows from calculations, the propagation of energetic nuclei through a TEPC is 
accompanied by production of fast (^-electrons which may escape the TEPC sensitive 
volume. This effect reduces the energy deposited to the TEPC, but it becomes less 
important for the same TEPC impacted by secondary nucleons. 

• Proper modeling of nuclear fragmentation reactions is crucial for describing 
microdosimetric distributions both on the beam axis and far from the beam. The 
nuclear fragmentation models of Geant4, G4QMD and G4BIC, are equally suitable 
for describing general features of the spectra, but both underestimate the fluxes of 
protons and neutrons far from the beam. This indicates the necessity to improve 
nucleus-nucleus collision models in calculating the angular and energy distributions 
of secondary nucleons and nuclear fragments. 
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• The relation between y~f and LET is not straightforward for nuclei which traverse 
a TEPC placed in water. Beam nuclei participate in fragmentation reactions in 
surrounding media. Thus they are partially replaced by fragments with reduced 
Z 2 . This reduces y~f measured in the presence of water. 

• The values of y~f and ya for a TEPC under direct impact of projectile nuclei are 
sensitive to the beam profile. 

• The MCHIT model is able to describe the spatial distribution of the total dose 
in the water phantom despite of six orders of magnitude decrease of dose with 
increasing distance from the beam axis. 

• The MCHIT model describes well the yd(y) distributions in the PMMA phantom 
irradiated by quasi-monoenergetic neutrons. This justifies the use of this model for 
calculating the TEPC response to secondary neutrons from 12 C beam. 

• The contribution of secondary neutrons to the out-of- field dose from 12 C beams is 
estimated from MCHIT simulations. Since experimental identification of neutrons 
is only possible with quite bulky detectors to be placed outside the water phantom, 
such microdosimetry measurements remain the only solution to estimate the upper 
limits for the dose from neutrons and their radiation quality. 

Our results can be extended to other measurements with TEPC, in particular with 
other kinds of nuclear beams, as the TEPC model (Far West Technology Inc., LET-1/2) 
we simulated in this work is used worldwide. A detailed consideration of contributions 
of specific nuclear fragments to yd(y) distributions calculated with MCHIT will be given 
elsewhere (Burigo et al. 2012). 
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